Hãy tưởng tượng một hệ thống vật lý—một khoản nợ đang tăng dần, một vật thể rơi tự do, hoặc một quần thể loài nguy cấp. Cấu trúc của phương trình vi phân cấp một (ODE) là cầu nối toán học cho phép chúng ta dự đoán trạng thái tương lai của các hệ thống này. cấu trúc của phương trình vi phân cấp một (ODE) là cầu nối toán học cho phép chúng ta dự đoán trạng thái tương lai của các hệ thống này. Nó định nghĩa mối quan hệ giữa biến độc lập $t$, biến phụ thuộc $y$ và tốc độ thay đổi tức thời của nó.
1. Phân loại cấu trúc
Ở cốt lõi, một phương trình vi phân cấp một liên hệ đạo hàm với các biến số: $$\frac{dy}{dt} = f(t, y) \quad (1)$$ hoặc dưới dạng ẩn $F(t, y) = 0$. Các phương trình được phân loại theo "khung xương" của chúng:
- Cấu trúc tuyến tính: Các phương trình như $\frac{dy}{dt} = -ay + b$ (2), trong đó hàm số là tuyến tính theo $y$. Lưu ý: Do đó, chúng ta chỉ sử dụng thuật ngữ 'nghiệm tổng quát' khi nói về các phương trình tuyến tính.
- Cấu trúc tự chủ: Khi tốc độ phụ thuộc hoàn toàn vào trạng thái hiện tại, $dy/dt = f(y)$. Những phương trình này thường có một ngưỡng giới hạn (T): mức dân số quan trọng mà nếu thấp hơn mức này, loài không thể sinh sản và sẽ tuyệt chủng.
- Cấu trúc chính xác: Được kiểm tra bởi điều kiện $M_y(x, y) = N_x(x, y)$. Nếu điều kiện này không thỏa mãn, như trong Ví dụ 3, thì không tồn tại hàm tiềm năng $\psi(x, y)$ nào thỏa mãn hệ phương trình.
Bước 1: Xây dựng mô hình
Các tình huống vật lý, ví dụ như VÍ DỤ 4 | Tốc độ thoát (một vật thể khối lượng $m$ được ném từ Trái Đất), phải được chuyển đổi sang ngôn ngữ toán học. Chúng ta cần xem xét lực hấp dẫn và vận tốc ban đầu $v_0$.
Bước 2: Tính ổn định và sự tồn tại
Chúng ta dựa vào điều kiện Lipschitz: $|f(t, y_1) - f(t, y_2)| \le K|y_1 - y_2|$ để đảm bảo nghiệm tồn tại và duy nhất. Nếu thiếu điều kiện này, "cấu trúc" của bài toán có thể bị sai lệch hoặc có nhiều nghiệm.
2. Nghiệm và trực quan hóa
Mọi hàm khả vi $y = \phi(t)$ thỏa mãn phương trình với mọi $t$ trong một khoảng nào đó được gọi là nghiệm. Về mặt hình học, chúng ta biểu diễn nghiệm này như một đường tích phân. Đối với phương trình Bernoulli, chúng ta sử dụng phép thế $v = y^{1-n}$ để tuyến tính hóa cấu trúc.
🎯 Nhận xét quan trọng: Phương pháp Euler
Trong VÍ DỤ 1 (cân bằng khoản vay $S(t)$ với lãi suất 12%), các xấp xỉ rời rạc dùng phương pháp Euler $y_{n+1} = y_n + f_n \cdot (t_{n+1} - t_n)$ thường lớn hơn các giá trị liên tục thực tế. Điều này xảy ra vì đồ thị nghiệm là lồi xuống, khiến các xấp xỉ bằng đường tiếp tuyến nằm phía trên đồ thị.
$\frac{dS}{dt} = rS - k \implies y_n = \rho^n y_0$